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Abstract 

To determine the nature of dark energy from observational data, it is important 
that we use model-independent and optimal methods. We should probe dark en- 
ergy using its density (allowed to be a free function of cosmic time) instead of its 
equation of state. We should minimize gravitational lensing effect on supernovae 
by flux-averaging. We need to include complementary data (for example, from the 
Cosmic Microwave Background [CMB] and large scale structure [LSS]) in a con- 
sistent manner to help break the degeneracy between the dark energy density and 
the matter density fraction. We should push for ambitious future supernova surveys 
that can observe a large number of supernovae at the highest possible redshifts. I 
discuss these and other issues that will be important in our quest to unravel the 
mystery of the nature of dark energy. 

Current supernova, CMB, and LSS data already rule out dark energy models with 
dark energy densities that vary greatly with time; with the cosmological constant 
model providing an excellent fit to the data. A precise measurement of dark energy 
density as a free function of cosmic time will have a fundamental impact on particle 
physics and cosmology. 



1 Introduction 



The data of type la supernovae (SNe la) reveal that the expansion of our 
universe is accelerating [Riess et al. 1998,Perlmutter et al. 1999]. This implies 
the existence of dark energy. Solving the mystery of the nature of dark energy 
will have revolutionary implications in particle physics and cosmology. 

The first most fundamental question about dark energy is whether the dark 
energy density depends on time. The time dependence of the dark energy 
density can be used to distinguish the various models. [Freese 2004] 

A powerful probe of dark energy is SNe la, which can be calibrated to be good 
cosmological standard candles [Phillips 1993, Riess, Press, & Kirshner 1995, Branch 1998], 
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and used to measure how distance changes with redshift. For a flat universe, 
the luminosity distance is 



dL{z)=cH,\l + z)J^ r 



(z'Y 



e(z)= n m (i + zf + (i - n m ) Px (z)/ Px (o) 



1/2 



(i) 



where px(z) is the dark energy density. Note that dark nergy only appears 
in the form of the function E(z); px{z) is on the same footing as the matter 
density fraction tt m in the context of observational tests. The dark energy 
equation of state is [Wang & Garnavich 2001] 



pxjz) (1-Kz) Px( z ) , 
Wx\Z) = — — = — — - 1. 
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2 Probing Dark Energy Using Its Density Instead of Its Equation 
of State 



Due to the integrals and parameter degeneracies in the distance-redshift re- 
lations of standard candles (see Eq.[l]), it is very difficult to detect the time- 
dependence of the dark energy equation of state wx (z) in a model- independent 
manner from realistic data [Maor, Brustein, & Steinhardt 2001,Barger & Marfatia 2001]. 
By measuring the dark energy density px{z), instead of wx{z), we effectively 
remove one integral, and increase the likelihood to detect the evolution of 
dark energy. [Wang & Garnavich 2001,Tegmark 2002] This would allow us to 
determine whether dark energy is vacuum energy, or some other form of exotic 
energy. 

It is easier to extract px( z ) from the data than to extract wx(z). This is easy 
to see mathematically: to obtain px(z), take a single derivative of data (see 
Eq.[l]); to obtain w x (z), we need to take a second derivative (see Eq.[2]). 1 

Fig.l explicitly demonstrates the advantage of the px{z) parametrization 
over the wx(z) parametrization, using recent observational data. Fig.l shows 
px{z) and wx(z) measured from 192 SNe la from the Tonry/Barris sample 



Note that in our analysis, no actual derivatives of data are taken; the best fit 
model is found by using a likelihood analysis based on Markov Chain Monte Carlo 
(MCMC). However, the mathematical relation between px{z) or wx(z) and the 
observable d^z) (through the apparent peak luminosity of SNe la) explains why 
px(z) should be better constrained than wx(z). 
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[Tonry et al. 2003,Barris et al. 2004], flux-averaged with Az = 0.05, and com- 
bined with CMB (shift parameter R = 1.716±0.062 [Spergel et al. 2003]) and 
LSS data (growth parameter f = f(z = 0.15) = 0.51±0.11 [Hawkins et al. 2003,Verde et al. 2002]). 
The solid and dashed lines enclose the 68.3% and 95% confidence regions 
respectively. [Wang & Mukherjee 2004, Wang & Freese 2004] Note that px{z) 



w/z) 





Fig. 1. wx(z) and px( z ) measured from 192 SNe la from the Tonry/Barris sample, 
flux-averaged, and combined with CMB (shift parameter Rq = 1.716 ± 0.062) and 
LSS (growth parameter = f(z = 0.15) = 0.51 ± 0.11) data. 

is more tightly constrained than wx{z). [Daly & Djorgovski 2004] and [Huterer & Cooray 2004] 
found similar results using a different sample of SNe la (the Riess sample). 

While using px {z) in the likelihood analysis of data encompasses all possible 
dark energy models, using wx(z) implies px(z) > since px(0) > (see 
Eq. [2]). To allow for exotic origins of dark energy, we should not impose a 
positive prior on the dark energy density. This again argues in favor of using 
Px(z) instead of wx(z) to probe dark energy. [Wang & Tegmark 2004] 



3 Go Deep: the Optimal SN Observational Strategy for Probing 
Dark Energy Evolution 



To determine whether SNe la are good cosmological standard candles, we 
need to nail the systematic uncertainties (luminosity evolution, gravitational 
lensing, dust) [Aldering 2004]. This will require at least hundreds of SNe la 
at z > 1. This can be easily accomplished by doing a supernova pencil beam 
survey using a dedicated telescope, which can be used for other things simul- 
taneously (weak lensing, gamma ray burst afterglows, etc) [Wang 2000a]. 

The inclusion of a large number of SNe la at high redshifts is important for 
probing the time evolution of dark energy density. Fig. 2 shows that a shallow 
and wide survey (the survey strategy for SNAP as of 2001) compares poorly 
with a supernova pencil beam survey in constraining dark energy density at 
high z [Wang & Lovelace 2001]. 

Clearly, in order to probe time evolution of dark energy density, it is essential 
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Fig. 2. A shallow and wide survey (the original survey strategy for SNAP) compares 
poorly with a supernova pencil beam survey in constraining dark energy density at 
high z [Wang & Lovelace 2001]. 



to conduct a SN survey that is as deep as feasible, so that the largest possible 
number of SNe la at high redshifts can be observed. 2 



4 Data Analysis Method 



Here is a recipe for optimally extracting dark energy density px{z) from ob- 
servational data [Wang & Mukherjee 2004]: 

1) Obtain SN la peak brightnesses with uncertainties in flux. 3 

2) Parametrize px(z) as a free function (given by px( z i)/ Px(0), i = 1,2, N, 
and interpolated elsewhere). Additional parameters are fl m and H . 

3) Use flux-averaging statistics in computing x 2 of SN la data [Wang 2000b]. 

4) Add other data (CMB, LSS) to help constrain Q m .[Knop et al. 2003] 

5) Use Markov Chain Monte Carlo (MCMC) to find marginalized probability 
distributions of estimated parameters. 

I will highlight the need for flux- averaging of supernovae and the advantage 
of MCMC below. 



2 The current SNAP survey strategy has been revised to image all the SNe la 
upto a redshift of 1.7 [Tarle et al. 2002], similar to a supernova pencil beam survey 
[Wang 2000a,Wang & Lovelace 2001]. 

3 Flux averaging is only straightforward if the intrinsic distribution of 
SN la peak brightness is Gaussian in flux (consistent with current data) 
[Wang & Mukherjee 2004]. 
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4-1 Minimizing gravitational lensing effect on SNe la by flux-averaging 



Since our universe is inhomogeneous in matter distribution, the apparent 
brightness of SNe la can be altered by the intervening matter distribution. 
[Wang, Holz, & Munshi 2002] gives a universal probability distribution for 
weak lensing magnification of SNe la. Because lensing only redistributes the 
light (flux) from the SNe la, the effect of weak lensing on supernovae can be 
minimized by flux-averaging. Therefore, flux- averaging justifies the use of the 
distance-redshift relation for a smooth universe in the analysis of SN la data. 
[Wang 2000b]. 

For x 2 statistics using MCMC or a grid of parameters, here are the steps in 
flux-averaging: 

(f) Convert the distance modulus of SNe la, fJ>o(z), into "fluxes" 

F( Zj ) = 10-^fe)- 25 )/ 2 - 5 = f^p) 2 • (3) 



(2) For a given set of cosmological parameters {s}, obtain "absolute luminosi- 
ties", {L(zj)}, by removing the redshift dependence of the "fluxes", i.e., 

L{z j ) = di{z j \B)F{z j ). (4) 



(3) Flux-average the "absolute luminosities" {L^} in each redshift bin i to 
obtain 

1 N i N 



(4) Place L % at the mean redshift Zi of the i-th redshift bin, now the binned 
flux is 

F(z l )=T/d 2 L (z l \s). (6) 



The 1-cr error on each binned data point F*, erf, is taken to be the root mean 
square of the 1-cr errors on the unbinned data points in the i-th redshift bin, 
{F;} (j = 1,2,..., iV), multiplied by 1/VN (see [Wang 2000a]). 
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(5) Apply "flux statistics" to the flux-averaged data, I F(zi 




(7) 



where F^Si) = (d L (z\s)/Mpc) 2 . 

A Fortran code that uses flux-averaging statistics to compute the likelihood of 
an arbitrary dark energy model (given the SN la data from [Riess et al. 2004]) 

can be found at http : / / www .nhn.ou.edu/ ~ wang / SNcode/ [Wang 2000b, Wang & Mukherjee 2004, V 

4-2 Likelihood analysis using Markov Chain Monte Carlo (MCMC) 

In a likelihood analysis, we need to find the high confidence region of our 
parameter space given the data. Traditionally, this is done by using a grid of 
parameter values that span the allowed parameter space. 

Markov Chain Monte Carlo (MCMC) uses the Metropolis algorithm: 

(1) Choose a candidate set of parameters, s*, at random from a proposal 
distribution. 

(2) Accept the candidate set of parameters with probability A(s, s*); otherwise, 
reject it. 

For Gaussian distributed observables, the acceptance function 



For sufficient sampling, the true probability density functions (pdf's) are re- 
covered. See [Neil 1993] for a review of MCMC. 

MCMC is much faster than the grid method, and scales linearly with the 
number of parameters. It gives smooth probability density functions (pdf's) for 
estimated parameters since they receive contribution from all MCMC samples 
(each with a random set of parameter values). Public software containing 
MCMC is available [Lewis & Bridle 2002]. 



5 Status and Outlook 



Using the latest SN la data, the Riess sample [Riess et al. 2004], which con- 
tains 6 SNe la at z > 1.25, together with CMB and LSS data, [Wang & Tegmark 2004] 
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obtained the most accurate measurements to date of the dark energy density 
Px as a free function of cosmic time (see left panel of Fig. 3). The solid lines 
indicate the 68% (shaded) and 95% confidence regions of px( z ) parametrized 
by its values at z = 0.467, 0.933, 1.4 (splined elsewhere for z < 1.4). The thick 
and thin dashed lines indicate the estimated px{z) parametrized by its values 
at z — 0.7, 1.4 (splined elsewhere for z < 1.4). A powerlaw form of px{z) is 
assumed at z > 1.4 (where there are only two observed SNe la); this allows 
the CMB and LSS constraints to be included in a self-consistent manner. This 
is important since it is incorrect (although common practice) to include priors 
on fl m derived assuming px(z)/px(0) = 1- [Wang & Tegmark 2004] Further, 
px (z) derived assuming fixed values of Q m may be misleading because of the 
degeneracy between px( z ) an d £l m . All our results have been marginalized 
over fl m and H . 




Fig. 3. Dark energy density px as a free function of cosmic time estimated 
from supernova, CMB (shift parameter i?o = 1.716 ± 0.062), and LSS (growth 
parameter /q = f(z = 0.15) = 0.51 ± 0.11) data. Left panel: The latest 
SN data, with px parametrized by 4 and 3 parameter splines; with estimated 
tt m = 0.278 ± 0.061 and n m = 0.279 ± 0.053 respectively. Right panel: The lat- 
est SN data (Riess sample) compared to previous data (Tonry/Barris sample, with 
estimated tt m = 0.332 ± 0.062). 

We compare the px( z ) estimated from the latest SN data (Riess sample) to 
that from the previous data (Tonry/Barris sample) in the right panel of Fig.3. 
Note that the Riess sample gives a substantially tighter constraint on px( z ) 
than the Tonry/Barris sample, and somewhat smaller Q m . Although the two 
estimates overlap at the 68% confidence level, the Riess sample clearly favors 
a dark energy density less dependent on time (with the cosmological constant 
model giving an excellent fit to the data). 

It is interesting that current data already rule out dark energy models with 
px{z) that vary a lot with time. Future deep supernova surveys on a dedicated 
telescope will greatly increase the number of observed supernovae [Wang 2000a] 



7 



This will allow improved calibration of SNe la as cosmological candles, leading 
to precise measurement of the dark energy density as function of cosmic time, 
which will have historical impact on particle physics and cosmology. 
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